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Quasicrystals can be considered, from the point of view of their electronic properties, as being 
intermediate between metals and insulators. For example, experiments show that quasicrystalline 
alloys such as AlCuFe or AlPdMn have conductivities far smaller than those of the metals that these 
alloys are composed from. Wave functions in a quasicrystal are typically intermediate in character 
between the extended states of a crystal and the exponentially localized states in the insulating phase, 
and this is also reflected in the energy spectrum and the density of states. In the theoretical studies 
we consider in this review, the quasicrystals are described by a pure hopping tight binding model 
on simple tilings. We focus on spectral properties, which we compare with those of other complex 
systems, in particular, the Anderson model of a disordered metal. We discuss "strong" and "weak" 
quasicrystals, which are described by different universal laws. We find similarities and universal 
behavior, but also significant differences between quasiperiodic models and models with disorder. 
Like weakly disordered metals, the quasicrystal can be described by the universal level statistics that 
can be derived from random matrix theory. These level statistics are only one aspect of the energy 
spectrum, whose very large fluctuations can also be described by a level spacing distribution that 
is log-normal. An analysis of spectral rigidity shows that electrons diffuse with a bigger exponent 
(super-diffusion) than in a disordered metal. Adding disorder attenuates the singular properties 
of the perfect quasicrystal, and leads to improved transport. Spectral properties are also used in 
computing conductances of such systems, and to attempt to resolve the experimental enigmas such 
as whether quasicrystals are intrinsically conductors, and if so, how conductances depend on the 
structure. 

PACS numbers; PACS numbers: 71.23.Ft, 71.10 Fd, 73.20 At 



I. INTRODUCTION 

Inside a perfect quasicrystal, such as the Penrose tiling 
[l[ or the octagonal tiling shown in FiglU a freely wander- 
ing observer would find a small reoccurring set of local 
environments, as in crystals. The sequence of environ- 
ments, however, does not repeat periodically, in the case 
of the quasicrystal. The repetitivity property states that 
any given pattern is guaranteed to be repeated within a 
distance proportional to the linear size R of the pattern. 
Random structures are very different in this respect, in 
that the number of patterns increases rapidly with the 
size of the pattern, and in this case identical regions of 
size R will be typically separated by distances that are 
exponentially large for large R. This geometric property 
has implications for the spatial extent of wavefunctions, 
as we will see later. Since the repetition is not periodic, 
one does not have a Bloch theorem in the quasicrystal. 
That the quasiperiodic structure is not random is per- 
haps best seen by calculating its structure factor, which 
turns out to have sharp Bragg peaks, as in a crystal [1]. 
The structure factor of a quasicrystal can also exhibit 5- 
fold 8-fold or 10-fold symmetries that are forbidden for 
crystals. 

A wave packet representing an electron that enters 
such a quasiperiodic medium will be scattered by the 
quasiperiodic potential, giving rise to a complex super- 
position of wave vectors and corresponding phase shifts. 
What does a typical wavefunction look like, and what 
are the allowed energy levels ? Would the electrons con- 



duct an electrical current, and if so, what factors deter- 
mine the value of the conductivity in a real quasicrys- 
tal ? Such questions have been asked, since the dis- 
covery by Schechtman et al Q of real quasicrystals in 
1984. Some answers have been provided, by direct ex- 
perimental measurement, and in numerical calculations. 
However, many basic questions concerning the proper- 
ties of even the simplest two-dimensional quasiperiodic 
structures remain unanswered. 




FIG. 1: Finite piece of a quasiperiodic structure 



In this paper we present a review of some electronic 
properties of a single electron in a quasiperiodic solid in 
two dimensions. This is the lowest dimension in which 
new physics due to the interplay between order and dis- 
order is expected to arise. A well known example is that 
of disordered crystalline solids described by the Ander- 
son model where any disorder, no matter how weak, 
"wins out" over the periodic order and exponentially lo- 
calizes wavefunctions, in one dimension. In two dimen- 
sions, this tendency is less marked, and one has the no- 
tion of weakly and strongly localized wavefunctions. A 
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transition between localized-delocalized transition only 
occurs, in the Anderson model, in three dimensions and 
above. A similar situation could be expected to hold for 
the quasiperiodic Hamiltonians. The case of one dimen- 
sion is the most tractable, in that analytical methods 
can be combined with numerical solutions in order to get 
a good appreciation of the expected physics. However, 
once again, one dimension is special in that the quasiperi- 
odic potential is almost always a relevant variable (in the 
sense of the renormalization group). 

One of the principal symmetries of the quasicrystal is 
its invariance under a scale change - inflations and defla- 
tions (to be described in Sec. II). The scale invariance of 
quasiperiodic structures, and in consequence, the poten- 
tial seen by the electrons, must play an important role 
in determining the spatial dependence and the energy 
of the eigenstates. The characteristic singular features in 
the density of states and other electronic properties are in 
fact a result of this symmetry. There have been attempts 
to exploit the heirarchical symmetry of the tiling in order 
to obtain renormalization equations 0] for tight-binding 
models, but the method breaks down for the pure hop- 
ping models that we discuss here. In a different context, 
another RG scheme was introduced for a spin model on 
this tiling, in order to calculate ground state energy and 
local magnetizations Q. It is therefore necessary, in the 
absence of analytical tools, to take a numerical approach 
to the electronic problem. 

Tilings can be deterministic ("perfect") or random. It 
is of interest to study the effects of including geometric 
disorder since real quasicrystals may well be intrinsically 
disordered in this sense. Random tiling models have been 
introduced (see review in Q) as the basic templates for 
naturally occurring quasicrystals. In contrast to the per- 
fect tiling, the random tiling has a lower free energy due 
to the associated entropy, explaining why materials might 
choose to condense with this type of long range order. 
Careful X-ray diffraction experiments must be done Q 
to distinguish the random tiling from the perfect tiling, as 
Bragg peaks occur in the same positions - this shows that 
the geometrical constraints are locally strong enough to 
ensure that even in random tilings the structure retains 
long range order with an infinite correlation length. 

The experimental motivation for our studies comes 
from a large number of intriguing results for the electrical 
conductivity of quasicrystals. Firstly, the low temper- 
ature electrical conductivity of quasicrystals is anoma- 
lously low, many orders of magnitude smaller than that 
of its constituent metals, aluminium, zinc, copper or iron, 
etc. Secondly, the conductivity rises with temperature, 
contrarily to the usual metallic behavior. Thirdly, when 
the structural disorder initially present in as-quenched 
samples of AlFeCu is removed by annealing, the conduc- 
tivity gets smaller . Similarly, in samples of AlPdMn 
prepared by different methods, the sample of hi ghe r 
structural perfection has the lower conductivity ll|, Il2l | . 



Disorder appears to facilitate transport in the quasicrys- 
tal, contrarily to its effect in disordered metals. Finally, 
studies of complex alloy system for the AlPdMn with very 
large unit cells show that the conductivity was higher for 
these than for a single grain quasicrystal of similar com- 
position, although both alloys had comparable structural 
quality. This could indicate that the transport is not just 
determined by local properties, but also by the long range 
quasiperiodic structural order. 

Many phenomenological Ts'l and numerical 14 1 studies 
have therefore addressed the problem of explaining the 
magnitude and temperature dependence of the electrical 
conductivity in realistic models. However studies of even 
the simplest of tight binding models show that the prop- 
erties of the quasicrystal are far more complex than those 
of crystals. Therefore, our focus in this paper is restricted 
to the simplest of two dimensional quasiperiodic tilings 
and determining for these whether states are expected to 
be localized or extended, whether electrons will conduct, 
and how best to describe the complex dynamics in such 
systems. 



I. DEFINITIONS AND GENERAL 
BACKGROUND 

We consider tight-binding Hamiltonians of the form 

in terms of onsite electron creation and annihilation op- 
erators c\ and Ci on a variety of two dimensional struc- 
tures. The first term allows for hopping between con- 
nected sites with amplitudes {ti\, while the second term 
allows for variations in the onsite potential. Although 
this Hamiltonian is grossly oversimplified, it contains the 
essential information about the quasiperiodic structure, 
namely, its geometry. For a discussion of the tight- 
binding Hamiltonian and its properties, particularly for 
one-dimensional quasicrystals, the reader is referred to 
[lit . In two dimensions, there have been many calcula- 
tions for this type of model and its dual version on the 
octagonal (Ammann-Beenker) and the Penrose tilings 

3 iH El, 111, 0, [m, m [23, im s [h, [23, m . 

We will now simplify the Hamiltonian even further by 
setting henceforth all of the hopping amplitudes equal 
to a constant, t = 1, and all the Vi ~ Q. This is be- 
cause we wish to focus exclusively on the effects of the 
quasiperiodic connectivity between atoms. Any varia- 
tions in the hopping amplitude, or of the local potentials 
has the effect of accentuating the quasiperiodic modu- 
lation and opening gaps in the spectrum. For example, 
if one includes onsite potentials V(Ti) that are propor- 
tional to the site coordination number, one gets, even 
for moderate values of V , bands separated by true gaps. 
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and wavefunctions which have a quasi one-dimensional 



character 23, 2 



The Anderson model for disordered conductors belongs 
in this class of tight-binding Hamiltonians. The disorder 
is realized by taking random hopping amplitudes and/or 
random onsite potentials or both. Such models are used, 
for example, to study critical properties at the metal- 
insulator transition in three dimensions. Spectral statis- 
tics in the Anderson model have also been discussed in- 
tensively, and compared with those of other classically 
chaotic systems. These studies ( see for example the re- 
view in [30| ) , along with the predictions of random matrix 
theory provide the background for the present discussion 
of the quasicrystal. 



Obtaining planar quasiperiodic structures 

A standard method of obtaining quasicrystals and 
their approximants consists of projecting down from a 
higher dimensional cubic lattice A subset of points 
belonging to an infinite strip of this hypercubic lattice 
is p rojected onto the physical plane. Figl2] (taken from 
[3l|) shows the projection onto the plane of a portion of 
the simple cubic lattice when the projection is rational 
(left-hand figure), and for an irrational orientation (right 
figure) where the structure never repeats. In these fig- 
ures, the different faces of the hypercube are projected 
onto tiles of different colors. A similar approach can be 
used to obtain the Penrose (projection from five dimen- 
sions), the octagonal tiling (projection from four dimen- 
sions) and the weakly quasiperiodic tiling described in 
Sec. Ill (projection of the simple cubic lattice). 

In calculations, the structures considered are often fi- 
nite periodic samples, called periodic approximants, ob- 
tained by rational cuts in the hyperspace. The unit 
cell of N sites can be made arbitrarily large, by modi- 
fying the plane of projection. The use of approximants 
has the advantage of eliminating surfaces and thus spu- 
rious energy levels - at the price of introducing a few 
defects in the structure. However, these are of dimin- 
ishing importance as the approximant size is increased, 
and the true quasicrystal structure is approached. For 
the square approximants of the octagonal tiling, we as- 
sume boundary conditions ^{x + L,y) — e'^''^^tl){x,y)] 



'ip{x,y + L) = e y 'ij}{x,y) where = 27m^/L and L is 
the repeat length of the approximant in each direction. 
An arbitrary external magnetic flux through a gauge field 
0, can also be introduced by suitably modifying the hop- 
ping amplitudes by a phase factor. As we will see, the 
addition of magnetic flux changes the symmetry proper- 
ties of the Hamiltonian, with a corresponding change in 
spectral characteristics. 

The N X N real symmetric matrix representing H is 
numerically diagonalized for each {kx-,ky} and the set of 
eigenvalues is used to compute the density of states and 





FIG. 2: a)a periodic projection and b) an aperiodic projection 
of the cubic lattice 



interesting statistical properties of energy levels. For pur- 
poses of computing the density of states, assuming large 
enough systems, it is in fact sufficient to compute the en- 
ergy spectrum at fc = 0. The two main statistical quan- 
tities considered are the nearest neighbor level spacing 
distribution P{s) and the level rigidity function S]^(i?). 
These quantities are introduced now, along with some 
known results and their interpretation. 



Level spacing statistics 

After ordering the energy levels corresponding to a 
given k and parity, we can compute the set of spacings 
Sn = En+1 — En- The mean level spacing is denoted by 
A. 

If the density of states p{E) — N^^ J2n ^{E~En) were 
constant as a function of the energy, one would simply 
go on to construct the histogram representing the prob- 
ability distribution P{s) defined by Eq[31 However, in 
general, the density of states does depend on the energy. 
In this case, the spacings must be redefined or "unfolded" 
so as to eliminate this trivial dependence on the position 
of the levels in the band. A systematic way to redefine 
spacings consists of fitting the integrated density of states 
by a smoothed "classical" part. This smoothed curve is 
then used to define a set of unfolded spacings (see [s^]). 
P{s) is now determined from the set of unfolded spacings 
{s„} using 
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TV 



(2) 



The issue of the statistical analysis of energy levels was 
raised first for the spectra of large nuclei. The first the- 
oretical calculations of level spacing distribution in com- 
plex systems were carried out by Wigner and by Dyson, 
for random matrices having independently distributed 
gaussian random matrix elements. The result depends 
on the class of matrix [s^. Three cases were distin- 
guished: real symmetric matrices (gaussian orthogonal 
ensemble or GOE), complex unitary matrices (gaussian 
unitary ensemble or GUE) and, finally, symplectic ran- 
dom matrices (GSE) 
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Pwd{s) = Aps^e-"^^'" (3) 

where /3 = 1,2,4 for GOE , GUE and GSE respec- 
tively. Values for the constants are, in the GOE case, 
Ai = 7r/2 and Ci = 7r/4. These distributions are writ- 
ten for spacings that have been normalized so that the 
mean level spacing A = 1. An impressive array of com- 
plex systems turn out to conform to one of these three 
expressions of Pwd{s)'- heavy nuclei, chaotic billiards, 
correlated fermions in regular crystals, etc. The Wigner- 
Dyson distributions thus appear to be a universal feature 
of classically non-integrable systems where the motion is 
ergodic. 

For the Anderson model, in the metallic regime many 
studies have confirmed that the level spacings obey GOE 
statistics. If time reversal symmetry is broken, by adding 
an external magnetic field, or by including a magnetic 
flux then GUE statistics are found. The third class of 
problems with /3 = 4 has also been realized in metals 
containing spin-orbit scatterers. 



sample do not overlap, and spatial correlations are neg- 
ligible, there is no bar to two different states having the 
same energy. This is the reason that a discrete level spec- 
trum corresponding to an insulator can be distinguished 
visually from that of of a weakly disordered metal. In 
FigO three different kinds of spectra are illustrated (en- 
ergies have been scaled so that there is the same mean 
spacing between levels for all three examples). In the in- 
sulator (left hand figure), distances between the energy 
levels of the insulator fluctuate greatly, and levels can ap- 
proach arbitrarily close to each other. The middle figure 
represents energy levels of a disordered metal and was 
obtained by diagonalizing the Anderson model Hamilto- 
nian. Here, levels repel each other due to non-negligible 
wavefunction overlap, and so there are fewer very small or 
very large spacings. Finally, the quasicrystal spectrum is 
intermediate between the first two examples (right hand 
figure) . 

Other probability distributions can be found in the 
literature, corresponding to other physically interesting 
problems. In cases where the classical dynamics of the 
electron is integrable, such as the motion of electrons in 
a crystal, the probability distribution of the set of all 
spacings is governed by a different law, the Poisson dis- 
tribution. 



Pp(.s) = e-^ (4) 

This behavior is found, as the name indicates, when- 
ever the levels are randomly distributed according to a 
Poisson probability (the left hand figure in Figl3]). It is 
the level spacing distribution in the strong localization 
regime of the Anderson model. At the critical point of 
the metal-insulator transition, there have been proposals 



3J] for another form of P{s) which is a hybrid between 
the exponential dependence at large s of Pp and the small 
s level repulsion of Pwd'- written above: 



(a) (b) (c) 



FIG. 3: Examples of spectra a) Uncorrelated energy levels 
b) disordered metal energy levels c) octagonal tiling energy 
levels 

A feature of note in the level spacing distributions 
Pwd is the fact that they tend to as s — > 0. This 
expresses the unlikelihood of levels being degenerate in 
this system - the phenomenon of so-called level repul- 
sion. Level repulsion occurs when electron wavefunctions 
are extended in the available volume, and when there is, 
in consequence, a non-negligible spatial overlap between 
states. In systems where wavefunctions are strongly lo- 
calized, on the contrary, states in different regions of the 



P,(s) = A0s'5exp(-as(2-'') (5) 

where /3 has the values 1 or 2 or 4 as mentioned earlier. 

The semi-Poisson distribution Psp is a special case of 
this type of law, with j3 = ^ — 1. The name semi-Poisson 
comes from the fact that this is the distribution for the 
distances between the midpoints between levels given a 
set of randomly chosen energy levels. There is level re- 
pulsion, but then the probability rises to a maximum 
and decays exponentially, as in the Poisson distribution. 
At the critical point, wavefunctions are expected to have 
multifractal character. This type of wavefunction leads 
to level repulsion, unlike states in the localized regime, 
but the effect is not as strong as that in the metallic 
regime. and Such multicritical wavefunctions have been 
linked to the level statistics at the critical point of the 
Anderson model 35] . This type of law has been observed 
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to hold in the Anderson model if one averages the distri- 
bution over all boundary conditions 3^, 37 , H, 3^ , at the 
critical point of the Harper model 4l| and in a disordered 
two particle Hubbard model f4^ . Another example with 
such statistics are the pseudointegrable billiards, studied 
by Bogomolny et al |42{ . 



The spectral rigidity function 

As we mentioned in the last section, the system of 
energy levels displays a certain amount of rigidity. One 
measure of this rigidity is given by the variance of the 
number of levels in an energy interval of width E: 



(6) 



where N{E) is the total number of levels within an 
interval of width E. The average is taken over all posi- 
tions of E in the spectrum. In random matrix theory, 
the rigidity has been calculated to be a relatively slowly 
increasing function of E, with 



^RMT 



{E) = -^H27tE/A) 



(7) 



In comparison, the rigidity of a completely random set 
of levels, where the P{s) is Poissonian, one has a linear 
dependence 



^%iE) = E 



(8) 



At the critical point, the spectral rigidities are also 
linear for small E, but with a smaller slope 
(x < 1), since correlations between levels begin to kick 
in, and lead to level repulsion. The semi-Poisson form 
will be described in more detail in the third section, for 
the case of the weakly quasiperiodic tiling. 

The spectral rigidity function can be obtained from the 
two-level correlation function K{E') defined by K{E') — 
{N{E)N{E + E'))/{N{E)y - 1. One can check that 
T.^{E) = 2 dE'{E~E')K{E'). Argaman et al es- 
tablished a connection between the dynamics of an elec- 
tron moving in the medium and the spectral form factor 
K. Adapting their calculations to a situation where an 
electron diffuses with an exponent a defined in terms of 
the root mean square distance d explored in a time t by 
d ~ t"' , then the probability of return to the origin scales 
as p{t) ~ t^^''-"' where d is the spatial dimension. It can 
be shown (see for example [ii^) that p{t) is proportional 
to K{t)/t where K{t) = ^/ K{uj)e-"^*duj. Transform- 
ing back to energy variables the power law in time for K 
becomes a power law in the energy for S^, 



Thus, if the rigidity function does indeed follow such 
a power law, one can from it deduce the exponent a de- 
scribing quantum diffusion in the medium. For disor- 
dered metals, one has normal diffusion, cr = ^, whereas 
as we will see, the motion in the quasicrystal is superdif- 
fusive and cr > ^. 

In numerical calculations on disordered metals one usu- 
ally finds different regimes of behavior of the spectral 
rigidity depending on the energy. The crossover energy 
scale between low and high energies is given by the Thou- 
less energy, Exh — fi/tL, where II is the time taken 
for the electron to diffuse to the boundary of the sam- 
ple. This time depends on the diffusion exponent via 
~ L^/*^. For small energies E << Eth, the behav- 
ior of is governed by the long time dynamics, which 
is completely ergodic since the electron has reached the 
boundary a large number of times. At such energies, 
the RMT logarithmic dependence is thus expected. For 
E » Eth, however, one expects to find a power law de- 
pendence, corresponding to the diffusive electron motion 
at times shorter than II- 



Wavefunctions and quantum diffusion 

The wavefunctions of the tight- binding Hamiltonian [T] 
are solutions of i/?/;'"^ = En'ij^'"-\ The set of amplitudes 
V'i"'' (* — where N is the number of sites) is numer- 
ically evaluated for each of the energy eigenvalues En- 
To determine whether a wavefunction is localized or ex- 
tended, one usually computes the inverse participation 
ratio (IPR) defined by 



N 



(n)|4 



(10) 



T.^{E) - E 



da 



(9) 



namely, the second moment of the probability density. 
For a localized state, the IPR remains constant as N in- 
creases, since such a state involves a finite fixed number 
of "participating " sites. For the extended metallic state, 
where a finite fraction of the total number of sites par- 
ticipate, the IPR decreases as the inverse of N. In qua- 
sicrystals, an intermediate situation obtains, with IPR 
cx where (3 is close to but less than 1 [45[ . 

The IPR is in fact a very limited probe of the real space 
structure of the wavefunction in the quasicrystal. A typ- 
ical wavefunction has an average behavior that decays 
slowly, but also has extremely large fluctuations from 
one site to the next as the figure shows. If one calcu- 
lates the exponents that describe the scaling of the mo- 
ments of i^iE) with N, one finds that the exponents are 
not just multiples of each other - this is the property 
of multifractality. A complete description of such "crit- 
ical" wavefunctions can only be given by specifying the 
complete set of moments of such a wavefunction. Such 
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wavefunctions have been studied in some detail for one 
dimensional systems like the Fibonacci chain. They have 
also been studied at the critical point of the Anderson 



model and the Harper model [41 



The diffusion of a wave packet on such tilings can be 
studied by computing the time autocorrelation function 
Ci{t), where 



(11) 



which gives the probability of the electron staying at 
the initial site i as a function of time. The mean square 
displacement for an electron that was initially located at 
the site i is 



(12) 



Again, one can determine exponents describing the 
long time behavior of these quantities, Ci(t) cx t^^ and 
r{t) (x t"'. It is this exponent a that has been related 
(see previous subsection) to the power law behavior of 
the spectral rigidity function S^(£^). Contrarily to the 
case of random systems where self-averaging leads to the 
same exponent regardless of which site is chosen as the 
origin, in the quasicrystal, a and S depend not only on 
the wave packet energy but also the site chosen as the 
origin. 



Transport 

Thouless [13] proposed a measure of intrinsic conduc- 
tance that is based on the sensitivity of each level to 
changes of the boundary condition. The energy of a lo- 
calized state is not affected by a change of boundary con- 
dition, whereas the energy of an extended state is. Thus 
if a level shifts as a function of k, the magnitude of the 
shift is an indication of the capacity for transport of that 
state. One can specify a "band width" for each band by 
Wn ~ En{0, 0) — £'„(7r, 0). The Thouless number is given 
by 



gth{En 



En+1 — En-1 



(13) 



Another quantity that has been used as a quantitative 
measure of the transport property of random systems is 
the set of energy level curvatures c„ = 9^i?(fc)/9fc^|fc_+o- 

A theoretical expression for the distribution of c (para- 
metric level statistics) in the case of weakly disordered 
metals has been obtained in RMT (see ^48j ) : 



where /3 = 1, 2 or 4. This law has been verified for disor- 
dered systems (ioj . 

In the case of the quasiperiodic tilings, both the 
Thoulesss number and the curvatures have huge fluctua- 
tions, as we will see below. The form of the distribution 
does not agree with the RMT form. 



II. A CASE OF STRONG QUASIPERIODICITY: 
THE OCTAGONAL TILING 



In this section we will take up the analysis of the octag- 
onal tiling with and without disorder. The word "strong" 
distinguishes this tiling from the weakly quasiperiodic 
tiling, the GRT, discussed in the next section. The oc- 
tagonal tiling has an eight-fold symmetric diffraction pat- 
tern, compared to another strongly quasiperiodic tiling 
- the Penrose tiling, which has ten-fold symmetry. This 
is not an exact rotational symmetry in real space as in 
ordinary crystals, instead, it reflects the fact that any fi- 
nite domain of the octagonal(Penrose) tiling, occurs with 
equal probability in each of eight (respective five) possi- 
ble orientations. 



Samples with and without disorder 

FigHK) shows a perfect (disorder-free) square approx- 
imant and the six types of vertices present in it. The 
infinite tiling and its approximants are composed from 
squares and 45° rhombuses. The infinite octagonal tiling 
can be obtained by projection from a higher dimensional 
cubic structure. For numerical studies, it is preferable 
to work with samples that do not have a boundary, the 
periodic approximants of the octagonal tiling (5Q] . 




FIG. 4: a) A square approximant of the perfect octagonal 
tiling b) The six local environments 



a/3 



(&^ + c2)(2+/3)/2 



(14) 



Some details about the properties of the octagonal 
tiling may help in understanding the structure better. 
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For more details, the reader is referred to 51, 153, \5 



The infinite tiling and its approximants has six local en- 
vironments as shown in figlH These are labelled A,B,...,F 
and correspond to coordination numbers z = 8, 7, 6, 3 
respectively. It is easy to show using the cut-and-project 
scheme, for example, that the relative frequencies of oc- 
currence of each type of vertex are 



/a = A /s = A ^; fc = 2A 
foi = A ^ = fD2] Ie = 2A '^■,fp = X ^ 




10 12 



FIG. 5: An approximant and its inflated version, showing the 
old and new tiles 



where A = 1 + \/2. This irrational is fixed by the 
condition of eight-fold symmetry of the quasicrystal and 
its equivalent on the Penrose tiling is the golden mean 
T — (\/5 + l)/2. An important symmetry of the qua- 
sicrystal is its self-similarity under a discrete scale trans- 
formation called inflation(deflation). Fig[5] shows how to 
redefine the tiles or to decimate sites, so as to obtain an 
infiated version of the tiling. The density of sites on the 
infinite tiling is reduced by a factor A^ after an inflation. 
The rational approximants transform into each other un- 
der successive operations of inflation or deflation. 

In [iBi, we have diagonalized the tight-binding Hamil- 
tonian of Eq[T] for increasing system sizes, Ng — 
239, 1393, 8119, using an extended Lanczsos routine that 
yields all distinct energy eigenvalues. 

Figinishows a sample obtained after one performs a cer- 
tain number of operations called "phason flips" wherein 
one permutes the positions of a square and two rhom- 
buses in the interior of a hexagonal shaped region. Iter- 
ating this operation at randomly chosen locations on the 
square approximants, one obtains a structure that has 
new environments that are not present in the original 
tiling. (Note: in the two dimensional case, contrarily to 
three dimensions, such randomizing leads to a reduction 
of the peak intensities in the structure factor [^). 
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FIG. 6: a)a square sample of the randomized octagonal tiling 
b) a tile flip (or phason) operation 



Density of states of the octagonal tiling 

Typical histograms representing the density of states 
(DOS) for the randomized and for the perfect octago- 
nal tiling are shown in the figures [7] and [H (In figlT] a 
delta function at E = corresponding to localized ring- 
like states around the eight-fold symmetric sites has been 
subtracted). Both graphs show a jagged dependence on 
the energy, with bigger fluctuations in the perfect case. 
The difference between the two DOS is not merely quan- 
titative, but qualitative. The fluctuations in the random- 
ized sample are bounded, and do not give rise to pseudo- 
gaps, where the DOS plunges to zero, as in the perfect 
system. When the DOS of samples of different sizes are 
plotted, it becomes evident that there is no smooth be- 
havior of the DOS in the perfect case: fluctuations are 
strong as N is varied at a fixed energy, and vice versa. 
This behavior is investigated in more detail below, where 
we consider the energy spacin gs fo r the two types of sys- 
tems. A multifractal analysis [26| gives a fractal dimen- 
sion of the DOS equal to one, for the perfect tiling, which 
thus has a single band. 



The level spacings distribution function P{s) of the 
octagonal tiling 

In analysing spacings between energy levels, one has 
first to classify the levels in groups according to their 
symmetry property. In the case of the approximants, we 
have one exact symmetry, namely that by reflection with 
respect to one of the diagonals of the square, x = y. 
There are other symmetries - the inflnite tiling has a bi- 
partite structure and so the spectrum is an even function 
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PiE) 




boundary conditions corresponding to zero and nonzero 
magnetic flux respectively (as explained further below). 
The data points are shown along with the Dyson- Wigner 
GOE probability distribution functions corresponding to 
f3 — I and 2, as expected. 



P(s) 



1 

OB 
OS 
7 
6 
5 
O.i 



FIG. 7: DOS of a perfect sample (from [5|]) 
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FIG. 8: DOS of a randomized sample (from [5; 



of the energy E. This is broken by the approximants, 
but can be restored by, for example, quadrupling the unit 
cell so that the boundary conditions do not frustrate the 
wavef unctions. In addition, the approximants have an 
approximate rotational symmetry [45|. In the calcula- 
tions mentioned below, levels were grouped according to 
even(odd) parity under reflection, but not separated into 
any further subgroups. 

In most situations, it is possible to determine the un- 
folded spacings without any ambiguity by separating the 
underlying smooth energy dependence of the DOS. In the 
case of the randomized tiling, for example, we are able 
to define a locally smooth DOS and use this to renor- 
malize spacings. The distribution obtained in 26] for 
these renormalized spacings is shown in FiglHl for two 



FIG. 9: P(s) for the randomized octagonal tiling with and 
without inclusion of a magnetic flux (from [sj]) 



We found, furthermore, that the same GOE distribu- 
tion is obtained when one introduces energetic (diagonal) 
disorder on the perfect octagonal tiling a la Anderson. 

Moving on to the case of the perfect tiling, the density 
of states fluctuations occur on all energy scales, so that 
it is not obvious how to separate out a "classical" part. 
The same type of problem occurs at the critical point of 
the Anderson and the Harper models. There, it is pos- 
sible to perform an averaging over all possible boundary 
conditions, and determine a local value of the DOS. This 
was used to unfold the spacings and thus obtain the crit- 
ical level statistics. In the case of the quasicrystal, we 
first calculated the distribution of the bare spacings, i.e. 
without unfolding. The bare spacin gs t urn out to follow 
quite well a log-normal distribution [26| 



Pis) 



2ttsB 



-(lns-lnso)^/2S 



(15) 



where the peak position, Inso, shifts to the left as the 
size of the tiling is increased. This is, of course, expected 
as the mean level spacing decreases A « W/N where W 
is the band width. FigfTOl shows the data for three sys- 
tem sizes, after the distributions were shifted so that the 
peaks coincide. Note that, unexpectedly, the width of the 
distribution, B, does not depend on the system size. An 
explanation is given in the section below where we outline 
the relation between the bare spacing and the unfolded 
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spacing distributions. There are deviations from the log- 



normal curve: for s small the distribution found is linear, 




FIG. 10: -P(s) for the perfect tiling for bare spacings for 
A'^=239 (stars), 1393 (triangles) and 8119(diamonds) (curves 
were shifted to that the peak positions coincide) 



FIG. 11: Statistics of levels for three cases: no level-unfolding, 
partial unfolding and maximal unfolding (from [Hj) 



The log normal -P(s) obtained for the bare spacings 
corroborates a fact already noted, namely the density 
of states fluctuates enormously in the perfectly ordered 
system. Next we consider possible ways of unfolding the 
spacings. The standard unfolding technique, involves a 
smoothing of the integrated DOS as described in (32i |. 
This smoothing has to be carried out on an energy scale 
that is small compared to the density of states fluctua- 
tions but large compared to the average nearest neigh- 
bor spacing, a condition that is not possible to respect, 
strictly speaking, in a quasicrystal. In an alternative 
definition, we divided the spacings by the local average 
value of the spacing defined over M consecutive levels: 
(s,) = {E,+M - E,_m)/^M. When M is large, then the 
renormalizations depend only weakly on the energy. In 
this case, the renormalized spacings obey a distribution 
close to that of the bare spacings. At the other extreme, 
for M = 1, the renormalization factor fluctuates strongly 
from one level to the next. In this case the renormalized 
spacings will obey a very different probability distribu- 
tion - this turns out to be the ubiquitous Wigner-Dyson 
distribution [s^l- As figure fT2l shows, one goes contin- 
uously from the log normal distribution to the Wigner- 
Dyson form simply upon changing the value of M in the 
unfol ding procedure. This was also confirmed by Zhong 
et al [27| using the standard unfolding procedure given 
in [s^l ■ A calculation of the level spacing distribution for 
a large but finite patch of the octagonal tiling confirmed 
that the same distribution is obtained in that case as 
well, despite the open boundary conditions [ist . 

The relation between the two types of level spacing dis- 
tribution is explained in terms of the heirarchical struc- 
ture of the quasicrystal in the next subsection. 



Relation between folded and unfolded statistics 

As the previous section showed, there does not appear 
to be a single way to define unfolded spacings for the 
infinite tiling, as the notion of "smoothing" is not well 
defined for this situation. In we have proposed a way 
to overcome this difficulty, which also resolves the issue 
of the relation between the two statistical distributions 
P{s) and P{s). When considering level statistics of ap- 
proximants, one can consider using the DOS of a smaller 
approximant to renormalize spacings of the next approx- 
imant in the series. Recall that approximants are related 
to each other by inflation/defiation transformations. In 
this way one can write the probability distribution of the 
bigger system in terms of a convolution involving the dis- 
tribution function of the smaller approximant. Iterating 
the process, and taking the initial distribution to be of 
Wigner-Dyson form, one can show that after a few iter- 
ations, the distribution for the spacings approaches the 
log normal form. In 55| we have analytically calculated 
the variance of the log normal distribution obtained in 
the large size limit, and we showed in particular, that it 
is independent of the system size, in agreement with the 
numerical calculations described earlier. 



The effect of adding magnetic flux 

The time reversal symmetry of the Hamiltonian in eq[T] 
can be broken by adding an external field, or, more sim- 
ply, by imposing a fiux via boundary conditions. This 
can be done by identifying two of the edges of the ap- 
proximant (similarly to when one wraps a graphene sheet 
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to form a nanotube) and considering a flux tube along 
the axis of the resulting cylinder. A magnetic flux (p 
along the axis of the cylinder will then create a total 
phase shift of e*^"^ on a closed loop around the flux 
tube. This phase shift is easily implemented by taking 
V'(x + L,y) = e"'^ip{x,y) i/^ix^y + L) = ip{x,y). The 
resulting distribution of spacings for the randomized sys- 
tem is GUE for sufficiently high values of (j>. This dis- 
tribution corresponds to a stronger suppression of small 
spacings than in the GOE case, since now P{s) ~ for 
small s (Eq[3]) . Fig. [5] shows the level statistics with and 
without magnetic flux, along with the theoretical GOE 
and GUE curves. 

In the case of the perfect tiling, one cannot apply this 
technique to create a magnetic flux through the samples 
because of the exact reflection symmetry of our samples. 
This symmetry causes the Hamiltonian to remain in the 
GOE class, because it is possible to find a real symmet- 
ric representation of H even after the introduction of a 
fiux tube through the sample ( i.e. if R is the operator 
for reflections (x, ?/ — s- y, x) then H is invariant under the 
combined operation RT where T is the time reversal op- 
erator). GUE statistics on the octagonal tiling can still 
be obtained by including an external magnetic field per- 
pendicular to the plane of the tiling in the Hamiltonian. 



The spectral rigidity function of the octagonal tiling 

Since level repulsion is present, wavefunctions are ex- 
pected to be extended in the quasiperiodic tiling. One 
can expect that a wave packet initially localized in some 
region will spread out over time. This is now shown by an 
analysis of the spectral rigidity. In the randomized tiling, 
the spectral rigidity function crosses over from the small 
E logarithmic law given in RMT to be T.'^iE) - \n{E) to 
a power law behavior. It is found that 



J:^{E) cx E'^ (16) 

in a large range of energies, with 7 w 1.7. The ex- 
ponent for quantum diffusion is therefore determined ac- 
cording to eqlHto be (T « 0.85. The perfect tiling shows a 
smaller value of the exponent, a w 0.8. The result for the 
perfect tiling is thus in good accord with Passaro et al 
|24l| who found an average value 0.78 for an electron dif- 
fusing in the octagonal tiling. Our value for the random 
tiling is higher than the value given by those authors of 
0.81. Both calculations conclude that, interestingly, dif- 
fusivity is higher on the disordered tihng. Disorder thus 
plays a role opposite to the one it plays in a crystal, in 
rendering the electron more mobile. 




FIG. 12: Spectral rigidity of the randomized tiling (solid cir- 
cles: GOE; open circles: GUE). The curves correspond to 
RMT predictions. Inset: log-log plot of the power law region. 
(from|2i]) 



Wavefunctions in the octagonal tiling 

Fig[T3] shows the probability distribution of an elec- 
tronic state near the band edge, for one of the square 
approximants. One can see the typical spiky spatial 
dependence, with some well-defined peaks, characteris- 
tic of critical wavefunctions. FiglH] shows results for 
p — P (participation ratio)/N(number of sites), plotted 
against the energy. Wavefunctions were investigated in 
detail on the Penrose tiling and are reviewed in [Zsj . 
where a similar plot is shown for the Penrose tiling. The 
data, in both the tilings, indicate that sites are more lo- 
calized in the center of the band, compared to the edges. 




FIG. 13: Probability densities on sites of an approximant of 
the octagonal tiling for a state of energy E m 3.98 



Transport in the octagonal tiling 

The results of calculations of the Thouless number in 
the perfect and the disordered octagonal tilings are shown 
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III. A FAMILY OF WEAKLY QUASIPERIODIC 
TILINGS 



FIG. 14: Plot of p =P(Participation ratio) /N for the octago- 
nal tiling for the approximant of 1393 sites 



in fig[T5] (taken from [5j]). The disordered tiling has, as 
expected, smaller fluctuations - note the change of scale 
between the figures. The perfect tiling appears to have 
the highest value of gth in the center of the band but 
the fluctuations here are also greater. A careful study 
of size dependence remains to be done for this system. 
Early calculations on approximants of the Penrose tiling 
[2^ showed that the Thouless number for that system in- 
creased towards the band edge. However, that model is 
dual to the vertex type model we consider here. Tsunet- 
sugu and Ueda also calculated the conductance of strips 
using the multichannel Landauer formalism, and found 
the spiky fluctuating dependence that is a consequence 
of the critical wave functions and spiky density of states. 
They found furthermore that the wavefunctions evolve 
from power-law-extended to power-law-localized, as en- 
ergy increased - since the tendency to open gaps increases 
with the energy in their model. Zijlstra has considered 
the octagonal tiling conductance by the same technique, 
but includes a potential energy term in the Hamiltonian 
which opens gaps in the spectrum. The results for the 
scaling of the conductance with size seem to indicate that 
despite the gaps, the wavefunctions for that system decay 
with power laws [2^ . 




FIG. 15: Thouless number for the perfect (left) and the ran- 
domized (right) tilings (from [s^ ]) 




FIG. 16: Rauzy tiling 
cubic lattice 



a two-dimensional projection of the 



Fig[TC] shows a structure introduced by Vidal and 
Mosseri called a generalized Rauzy tiling [56|(GRT). It 
is obtained by projecting a cubic lattice along an irra- 
tional direction. The sample shown is an approximant, 
i.e., it has been constructed so as to allow periodic bound- 
ary conditions, but the infinite structure is quasiperiodic. 
The tight binding hamiltonian for approximants is very 
simple to write down contrarily to the quasicrystals dis- 
cussed earlier. This is due to the fact that this lattice 
has a codimension of 1 (the codimension is the differ- 
ence D — d, where d is the physical dimension and D the 
dimension of the hypercubic lattice). 

For periodic approximants, the connectivity matrix is 
a band diagonal (Toeplitz type) matrix. The nonzero 
matrix elements are situated at distances of Fns, Fn-2 
and Fn-i from the diagonal, where the generalized Fi- 
bonacci numbers {Fn} are obtained from a three term 
recursion relation Fn = i^„_i +-F'„_2 + ^n-3, with the ini- 
tial conditions = 0; Fo = Fi = 1. The ratio F„/F„_i 
tends to the value a « 1.839..., solution of the equation 
= + X + 1. These matrices correspond to approxi- 
mants of increasing size as n is increased. For example, 
with n = 5 one has a 13 site periodic problem with a site 
connectivity matrix C given by 
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It is thus not only easy to write down the Hamihonian 
for big approximants for numerical calculations, but one 
has in this indexation of the sites a convenient basis for 
analytical calculations, as opposed to the octagonal and 
Penrose tilings. 

The question now arises: are there any physical con- 
sequences of reducing the codimension ? From the point 
of view of geometry it seems clear that the fluctuations 
of the geometry will be reduced when there are fewer 
"degrees of freedom" in the problem. As we saw in 
the preceding section, a study of the eigenvalues of the 
tight binding Hamiltonian can yield information about 
the wavefunctions and the dynamics of an electron in the 
medium. We have already seen the use of statistical tools 
with which to quantify the degree of complexity of the 



Hamiltonian. The results obtained for the GRT in [57 1 
gave a number of insights into the similarities and dif- 
ferences between this quasicrystal and the ones that had 
been studied before, as we will now describe. 

After diagonalization of the Hamiltonian, one obtains 
the spectrum and the DOS as shown in the figurefTTl The 
DOS rises sharply towards the band center, recalling the 
van Hove singularity of the square lattice (and one sees in 
figHSthat the GRT has many small patches of sites of co- 
ordination number 4). The DOS does not have the large 
fluctuations displayed by the octagonal tiling, instead, 
there appears to be an underlying smooth component to 
the GRT density of states. This enables us to compute 
the level statistics analysis using an unfolding technique 
similar to the one used for random systems. 



Level statistics of the GRT 

The distribution of spacings is not of the Wigner- 
Dyson type. We find instead a good fit to semi-Poisson 
statistics: 



FIG. 17: Density of states of the GRT 



Figs fTSl show the plot of the data obtained for the dis- 
tribution of level spacings, along with the theoretical semi 
Poisson and the Poisson curves. The agreement is good 
over a large range of s, although the data falls off faster 
than the theoretical prediction at large s. 

In the present case of the GRT, and the results for its 
P(s), one can speculate that the electron motion is not 
integrable as in a crystal, but not chaotic as in the oc- 
tagonal tiling. The GRT thus represents an intermediate 
case between the integrable square lattice and the non- 
integrable octagonal tiling. This conclusion is in keeping 
with our intuition that the quasiperiodicity is, so to say, 
weakened due to the small codimension. 

One expects, correspondingly, that the electron motion 
should approach the ballistic limit, with the root mean 
square distance travelled being linear in time t. This can 
be checked by computing the spectral rigidity function 
next. 



Spectral rigidity of the GRT 

The spectral rigidity corresponding to the semi Poisson 
case was calculated to be 



E2 



sP 



[£;+-(l-exp-4£;)] 



(19) 



Psp[s) = 4sexp 



-2s 



(18) 



which has a linear dependence at small energy, with 
a slope of one-half. Fig[Tn] shows the spectral rigidity 
(without unfolding) for the GRT. The behavior found 
agrees for small -E, where it is linear, with a cross over 
to quadratic. The latter behavior at large E (or short 
times) reflects the fact that the motion of the electron 
is close to ballistic, as we have already surmised and the 
exponent for quantum diffusion a k, \. 
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FIG. 18: P(s) of the GRT for three different sizes. Curves 
correspond to semi-Poisson and Poisson laws discussed in the 
text (from [57[) 



one goes out to the band edge. One can propose a hand- 
waving explanation of this unconventional behavior in 
the quasicrystal. In a disordered metal, the long wave- 
length (low energies, close to the band center) modes 
are relatively less scattered by the impurities, and so 
these wavefunctions remain long ranged and contribute 
to transport. For higher energies, the wavefunctions are 
more sensitive to scattering and will tend to become more 
localized as the wave vector (and the energy) increases. 
In the case of the quasicrystal, for very low energies, 
where the wavefunction modulation is slow, the lack of 
translational invariance is more strongly felt compared 
to higher energies, where the wavefunction varies on the 
scale of the small crystalline patches that we alluded to 
earlier. This leads one to expect that the quasiperiodic 
fluctuations of geometry are less effective in localizing 
the wavefunction at high energies than at low energies. 
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FIG. 19: Spectral rigidity of the GRT on log log scales (from 



An explicit calculation of the root mean square dis- 
tance have been carried out by Triozon et al [HI], for 
wave packets of varying energies. The effective exponent 
a found by these authors is about 0.9 at the center of 
the band, goes down to a minimum value of about 0.85 
and then increases again, approaching 1 as the energy 
increases towards the band edge. The linear behavior of 
at small energies agrees with the semi Poisson for- 
mula given in Eq[T9l The increased diffusivity at the 
band edges was found by Triozon et al [H^ to be accom- 
panied by an increase of the wavefunction participation 
ratio (PR) as energy increases. This is reminiscent of the 
behavior of the IPR on the Penrose tiling described in 
the previous section. 

One sees that the energy dependence of a behaves in 
the opposite way of what one expects in disordered sys- 
tems. There, wavefunctions tend to be most delocalized 
at the band center, and typically get more localized as 



IV. DISCUSSION AND CONCLUSIONS 

We have presented studies of the energy level statis- 
tics in several two dimensional quasiperiodic structures 
and compared them with periodic as well as disordered 
structures. 

The strong quasicrystals such as the octagonal tiling 
and also the Penrose tiling, are described by two different 
forms of the level spacing distribution. There is the dis- 
tribution P(s) of the maximally unfolded spacings, where 
fluctuations in the density of states have been compen- 
sated down to the smallest energy scale, has the univer- 
sal Dyson- Wigner form. These statistics are found irre- 
spective of whether a geometric disorder is added or not. 
The tilings resemble, in this respect, other complex sys- 
tems, with classically chaotic, ergodic trajectories, and 
the "strong" quasicrystal falls into this category. If one 
considers the bare unrenormalized spacings, -P(s) has a 
log normal form for both the octagonal and the Penrose 
[60| tilings, as in some other problems involving heirar- 
chical processes. In this case, the log normal spacing 
distribution reflects a heirarchical structure of the spec- 
tral density, which has huge fluctuations. The DOS for 
the randomized case, on the contrary, is closer to that of 
a typical disordered system. 

We then considered a weaker quasicrystal, the gener- 
alized Rauzy tiling, obtained by projection of the simple 
cubic lattice. Here the spacings have a different probabil- 
ity distribution - the semi-Poisson law. This tiling rep- 
resents an intermediate situation between the strongly 
chaotic and the integrable systems. The Rauzy tiling is 
in this respect closer to the crystal than the octagonal 
tiling, due to its smaller codimension. 

We showed that the dynamics of the electron, as de- 
duced from the spectral rigidity function also follows a 
different behavior in the strong and the weak quasicrys- 
tals: propagation with an average diffusion exponent of 
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about 0.8 in the perfect octagonal tiling, corresponding to 
a superdifTusive dynamics, while having a value close to 
unity for the generalized Rauzy tiling, corresponding to 
near-ballistic propagation. Disordering the perfect octag- 
onal tiling results in a slightly bigger value of the average 
exponent for anomalous diffusion. 

One can ask what occurs in the limit of increasing, 
instead of decreasing, the codimension of the two di- 
mensional tilings. This would result in structures hav- 
ing an increasing number of local environments, and in- 
creasingly stronger quasiperiodic disorder with D. The 
motion of the electrons is expected to become progres- 
sively more hampered. Destainvillc et al (soj have shown 
that there is no localization even in the limit of infinite 
D. As codimension is increased, the exponent for diffu- 
sion decreases, and appears to tend to the value a — 0.5 
i.e. ordinary diffusion. The level spacing distribution has 
not been calculated, but we would expect, based on the 
results for the other systems discussed here, this to be 
Wigner-Dyson. 

Finally, a word about the conductance properties of 
two-dimensional quasicrystals. The spectral studies de- 
scribed here lead us to expect that the pure hopping 
single component quasicrystal is close to being a metal- 
lic conductor. The quantum diffusion studies show that 
there are huge fluctuations as a function of the initial po- 
sition of the particles, their energy, and the time elapsed 
- leading to some effective average behavior that is su- 
perdiffusive. Similarly, huge fluctuations occur in the 
Thouless number for the perfect tiling, but are atten- 
uated in the randomized structures. The effect on the 
conducting properties is not well understood. This effect, 
as well as the detailed characteristics of this distribution 
remain to be investigated in detail. 
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the Orsay campus, for making these studies possible. 
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